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CN Abstract 

In this paper we explore a family of congruences over N* from which one builds a sequence of 
symmetric matrices related to the Mertens function. From the results of numerical experiments, 
we formulate a conjecture about the growth of the quadratic norm of these matrices, which implies 
the Riemann hypothesis. This suggests that matrix analysis methods may come to play a more 
important role in this classical and difficult problem. 
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^ 1. Introduction 

Among the numerous statements equivalent to the Riemann hypothesis, a few have been 
^ formulated as matrix problems. The Redheffer matrix A„ = (ai,j)i<ij<n is an n x n matrix 

T-H defined by a^j = 1 if j = 1 or if i divides j, and a^j = otherwise. R. Redheffer [16] proved that 

^ the Riemann hypothesis is true if and only if 

cn det{An) = 0(n^/^+'), for every e > 0. 

F. Roesler [17] defined the matrix Bn = {bij)2<i,j<n by bij = i — 1 if i divides j, and bij = — 1 
00 otherwise, and he thus proved that the Riemann hypothesis is true if and only if 

det{Bn) = 0(n!n~^/^+'), for every e > 0. 

Both the Redheffer and Roesler matrices are related, via their determinants, to the Mertens func- 
tion, which is by definition the summatory function of the Mobius fimction (sec [1] p. 91). These 
matrices are not symmetric, and many eigenvalues must be computed to estimate the determi- 
nants. 



The matrix that we introduce in this paper is also related to the Mertens function, and 
it is symmetric. In Section 2, we will prove that the Riemann hypothesis is true if 

||A^„|| = 0(ni/2+"), for every e > 0, 

where ||A^n|| denotes the quadratic norm of the matrix, i.e. its spectral radius. In Section 3, from 
the results of a series of numerical experiments, we will conjecture that ||A^n|| = 0{n^^'^~^^), for every e > 
0. Consequently, the search for a good estimate of ||A^n|| turns out to be a key problem. At this 
prospect, it should be noted that, contrary to both the Redheffer and Roesler matrices, only one 
eigenvalue must be estimated in the matrix Mn- 
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2. Construction of the matrices Adn 

2.1. A family of congruences over N* 

Definition 1. To each n G N*, we associate an equivalence relation Tl over N*, defined by 

i n j [n/i] = [n/j] . 

Example 2. For n = 16, TZ possesses the eight equivalence classes 

{1}, {2}, {3}, {4}, {5}, {6, 7, 8}, {9, 10, 11, 12, 13, 14, 15, 16} and {17, 18, • • • }. 

Due to the interval structure of these classes, we can unambiguously identify each class by 
its largest representative (with the convention that oo denotes the largest representative of the 
unbounded class). For clarity, the representatives are written in plain letters, and the classes are 
written in letters with a hat. We denote by S the set of the largest representatives and by S the 
set of the classes, i.e. S = W /TZ. We also set 5 = 5 \ {oo}. 



Throughout Section 2, most of the objects we define, such as the sets TZ,S defined above, 

and the matrix Ai introduced in Proposition 21, depend on the integer n. However, in order to 
simplify the notations in this section, we will not index these objects by n, since there is no risk 
of ambiguity. 

Example 3. For n = 16, 

5 = {1,2, 3, 4,5, 8,16, oo}, S = {1,2,3,4,5,8,16}, 5 = {1,2, 3, 4, 5,8, 16, oo}, 

T = {1}, 2 = {2}, 3 = {3}, 4 = {4}, 5 = {5}, 8 = {6, 7, 8}, 16 = {9, 10, 11, 12, 13, 14, 15, 16}, 

oo = {17,18, •••}. 

Proposition 4. Let n he fixed in W , and let S he defined as ahove, i.e. S is the set of the largest 
representatives of the classes ofTZ. For each k in S, we set k = [n/k]. 

1. For each k in S, we have k E S and k = k, which means that the map k ^ k is a decreasing 
involution on S. Actually, k ^ k is just the order reversing map on S. 

2. The set S can he descrihed precisely hy the following alternative : 

Ifn< [^/n\^ + [^/^ then S = |l,-- - ,[v^ = \^^... ,n = T}, 

hence = 2[y^] — 1; 
ifn > [^/n]'^ + [-v/n] then <S = |l, • • • , [^/n\, [^/n], ■■■ ,n = T| , 

hence = 2[^/n\. 

Proof. 1. Let k e S and k = [n/k], which means that kk < n < kk + k. That is to say, 

n , n 

^ <k< ^. 

k+l k 



Since A; is an integer, it follows that 



n 




n 




< 




k+ 1. 




J. 



which proves that k G S. 



Let k E S. Since 



n 






k + l 


< 




.k. 



, it follows that 



n 



k + l 



< 



k. That is to say, n < kk+k. 



n 



From this last inequality and the fact that kk < n, we deduce that k < = < k + 1, which 

k 

means that k = [n/k] = k. 
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2. We begin to prove that each singleton set {k}, with 1 < k < -^/n, is a class. Indeed, 

if k < \^/n\, then we have successively k + 1 < [v^lj i i = tji ^ > 

k k + 1 k[k + 1) 

This last inequality means that A; and A; + 1 do not belong to the same 

class, which is to say, k e S. 

Considering now the case k = [\/n\ , there are two possibilities : 
n n n 



n 






k+l_ 


< 




.k. 



(a) either 



(b) or 



k 
that 



k k + 1 k{k + l^ 
n n 

= —7- — < 1, hence 

fe+l k{k + l) 



> 1, hence 

n 



n 




"n" 


_k + l 


< 




-k- 



n 



< k < —, and since k is an integer, it follows 
k+1 ~ k ^ 



n 




n' 


_k + l_ 


< 




.k. 



In both cases, [-^/n] and [-^/n] + 1 do not belong to the same class, so [^/n\ G S. 
Now that we have proved that {1,2, • • • , [-v/n]} C S, let k £ S with k > [\/n]. Therefore, 
k satisfies the inequalities k > \fn^ n/k < \fn and k < [^/n\. In other words, [^/n\ is the 
largest element of S such that k < k. Using the fact that k kis a, decreasing involution on 
S, and distinguishing the two cases [■\/n\ = [■\/ri\ and [■\/ri\ < [■\/n\, we deduce the expected 
form of S. 

To conclude, we rewrite the condition [y^] = [^/n\. Set -^/n = k + a with k & N* and 
< a < 1. We have n = k'^ + 2ak + a^, i.e.: n/k = k + 2a + a'^/k, so the following 
equivalences hold : 

[^/n\ = [^/n\ 4^2a + a^/k <l-^ + 2ka - k<04^n<k'^ + k, 

and this completes the description of S. 

Remark 5. A synthetic formula for #5, one valid for all n G N*, is 

#5=[V^] + [V^^Tl74-l/2]. 

Lemma 6. For all positive integers n,i,j, we have 

= [n/ij] ■ 

Proof. Write n/i = u + a, with u £N and < a < 1. 

Hence [n/i]/j = u/j and n/ij = u/j + a/j, from which it follows that [[n/i] /j] < [n/ij]. If this 
last inequality were strict, then there would exist an integer v such that u/j < v < u/j + a/j. 
Hence, u<vj<u + a<u + l, which is impossible since both u and vj are integers. 



Proposition 7. TZ is compatible with the multiplication over W , meaning that, for all i,j, k G N*, 
we have i TZ j =^ ik TZ jk. Therefore, the formula ij = ij defines an induced multiplication over 
S (recall that i denotes the class of i, and oo the class of every integer strictly larger than n). 

Proof. Assume that i IZ j, so [n/i] = [n/ j]. Using the previous lemma we deduce 
[n/ik\ = ^n/i\/k] = [[n/j]/A;] = [n/jk\, which is to say, ik TZ jk. 
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The set N*, equipped with the usual multiphcation, is a commutative semigroup. Since TZ 
is compatible with the multiplication, the quotient set S = N* /TZ, equipped with the induced 
multiplication, is also a commutative semigroup. Prom now on, for each fixed n in N*, the 
equivalence relation TZ will be called a congruence, and any two integers i,j such that iTZj will be 
said to be congruent. 

Example 8. Multiplication table of <S = W/TZ, for n = 16. 





1 


2 


3 


4 


5 


8 


16 


oo 


1 


1 


2 


3 


4 


5 


8 


16 


oo 


2 


2 


4 


8 


8 


16 


16 


oo 


oo 


3 


3 


8 


16 


16 


16 


oo 


oo 


oo 


4 


4 


8 


16 


16 


oo 


oo 


oo 


oo 


5 


5 


16 


16 


oo 


oo 


oo 


oo 


oo 


8 


8 


16 


'55 


'55 


'55 


'55 


'55 


oo 




















16 


16 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 



2.2. Three Z-algebras 

The Z- algebra of a semigroup G is the set Z'^ of maps from G to Z equipped with the 
convolution product -k, defined naturally as follows : if a and b are elements of Z*^, then c = a-kh 
is the map defined by 

Mt eG, c{t)= a{r)b{s). 

r,s&G:rs=t 

Of course, this makes sense if the above sum is finite, a condition that is always satisfied in the 
remainder of this paper. 

2.2.1. The algebra A of the semigroup N* 

The algebra A = Z^* of the semigroup N* is the algebra of Dirichlet series (with integer 
coefficients) equipped with the convolution product also called Dirichlet product (see [1] p. 29). 
This algebra possesses some well-known properties : 

Properties 9. 1. If a = (oi, • • • , afc, • • • ) and b = (bi, ■ ■ ■ ,bk, • • •) are elements of A, and 
c = a*b, then = Ylij=k o-i^j- 1^ particular, for Cj, ej, the i-th and j-th vectors, respectively, 
of the canonical basis of A, then we have * ej = eij . 

2. The unit of A is ei = (1, 0, • • • , 0, • • • ). An element a = (ai, • • • , an, • • • ) is invertible if and 
only if ai = ±1. 

3. The inverse o/u = (1, 1, • • • ,!,•••) is fi, the Mobius sequence (see [1] p. 31). 

2.2.2. The algebra A of the semigroup S, the quotient algebra A. 

If we consider the semigroup S = N*/7^ equipped with the induced product, then its algebra 
A = Z"^ is a Z-algebra of dimension #5, for which 5 is a basis. However, more interesting for 
our purposes is the quotient algebra A. = A/ooA where 00^4 = Zoo is the principal ideal of A 
generated by oo. Let w be the canonical projection of A onto A.. We print in bold the images 
under w of the vectors of the basis 5, and more generally any vector in A. For instance, w{k) = k, 
ro(oo) =0. As /c runs through the set S (see Example 3), k runs through a set denoted by «S, 
which is a basis of A, called the canonical basis of A. Of course #«S = #5, a quantity that has 
been computed in Proposition 4. Using these notations, it is easy to construct the multiplication 
table of the basis S, from the multiplication table of 5, as follows : 
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1. Remove the last row and the last column from the table of S, and replace the remaining 
symbols oo by (this expresses the fact that ro(oo) = 0). 

2. Remove the hats and rewrite the integers in bold letters (this expresses the fact that w{k) = 
k). 

Example 10. Multiplication table of the canonical basis S, for n = 16. 





1 


2 


3 


4 


5 


8 


16 


1 


1 


2 


3 


4 


5 


8 


16 


2 


2 


4 


8 


8 


16 


16 





3 


3 


8 


16 


16 


16 








4 


4 


8 


16 


16 











5 


5 


16 


16 














8 


8 


16 

















16 


16 





















2.3. A natural morphism from A to A. 

Definition 11. We call ^9 the morphism of Z-modules defined by 



A^A 



, where denotes the i-th vector of the canonical basis of A. 



Proposition 12. ■& is a morphism of "L-algehras. That is, for all a, b in A, we have 

7?(a*6) = i?(a)i?(6). 

Proof. Since is linear, we have only to prove the result for arbitrary basis vectors Cj, ej. Using 
Propositions 7 and 9, we have '&{ei -k ej) = i!}{eij) = ij = ij = 'd{ei)'d{ej). 

Proposition 13. The map tt = w o 'd from A to A. is a morphism of "L-algehras. 
The image of a = (ai, • • • , a^, • ■ ■) & A by the morphism tt is 



Proof. The only thing to check is that the sum J^K^it ^^^^ defined. This results from the 
fact that, for every A; G <S, A; is a finite subset of N*. 

Corollary 14. The m,orphism, tt has the following properties : 

1. It maps M = (1, 1, • • • , 1, • • • ) to u = {#k)k£S. If for each k G S we denote by k~ the 
predecessor of k in S, with the convention that 1~ = 0, then ^k = k — k~ . Hence, 

u = {k- k~)k&s- 

2. It maps fx = {fx{l),n{2),- ■ ■ ,/x(fc),---) to /j, = ^E^efe M/^))^^^- Hence, 

tJ,= {M{k)-M{k-))^^^, 
where M denotes the Mertens function. 
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Example 15. For n = 16 wc have 

1. u = (1,1,1, la, 3,8) = 1 + 2 + 3 + 4 + 5 + 3 8 + 8 16. 

2. fi = (1,-1,-1,0,-1,0,1) = 1- 2- 3- 5 + 16 

= ^(1)1 + fi{2)2 + //(3)3 + /i(4)4 + ;u(5)5 + (M(8) - M(5)) 8 + (M(16) 

2.4- The regular representation of the algebra A. 
For every a G the map 



M(8))16. 



X I- 



ax 



is hncar, and it is represented in the canonical basis «S of ^ by a matrix />(a). We denote by 
s = #S the dimension of A {s was computed as a function of n in Proposition 4), and we denote 
by MsC^), the algebra of square matrices of size s with integer entries. The map 



P ■ 



a I 



is called the regular representation of A (see [7] p. 56). Moreover, this representation is faithful, 
i.e. the morphism p is injective. The set of all matrices p(a), for a G is therefore a commutative 
sub-algebra of Ais{'^), of dimension s, for which the matrices p(k), k G «5, form a basis. Finally, 
since there is a natural bijection between S and S, we choose S as the indexing set for the rows 
and the columns of the matrices p(a). For instance, when n = 16, the last column of any matrix 
p(a) does not have index 7, but index 16. 

Example 16. For n=16, in addition to p(l), which is the identity matrix, the matrices repre- 
senting 2, 3, 4, 5, 8, 16 (where most of the zero entries are left blank for legibility) are 



P(2) 

1 
2 
3 
4 
5 
8 
16 



2 3 4 5 8 16 



1 1 



P(3) 
1 
2 
3 
4 
5 
8 
16 



1 2 3 4 5 8 16 



1 1 1 



P(4) 
1 
2 
3 
4 
5 
8 
16 



P(5) 



1 
2 
3 
4 
5 
8 
16 



1 2 



1 

2 
3 
4 
5 
8 

16 



:i i 



5 cS 



/;(16) 



1 

2 
3 

4 
5 
8 
16 



1 2 3 4 5 8 16 



1 1 



1 2 



8 16 



Furthermore, following Example 15, the matrices representing u and fj, are 



P(u) 


1 


2 


3 


4 


5 


8 


16 


p(m) 


1 


2 


3 


4 


5 


8 16 


1 


1 














1 


1 












2 


1 


1 












2 


-1 


1 










3 


1 





1 










3 


-1 





1 








4 


1 


1 





1 








4 





-1 





1 






5 


1 











1 






5 


-1 











1 




8 


3 


2 


1 


1 





1 




8 





-1 


-1 


-1 





1 


16 


8 


4 


3 


2 


2 


1 


1 


16 


1 


-1 


-2 


-1 


-2 


-1 1 
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on which we verify that for a G the coefficients of a in the basis S appear in the first column 
of p(a), cf. Example 15. 



Definition 17. For n G N*, let T be the symmetric matrix of size s = whose entries are all 
I's above the second diagonal and O's strictly below. 



Example 18. For n = 16, the matrix T is 



111111 
11111 

1111 
1 1 1 
1 1 
1 



Lemma 19. For every k G <S and every i,j G <S, the following equivalence holds : 

{Tp{k))^. = 1 ^ ij < [n/k] ^k< [n/ij]. 

Proof. Let k G «S, j € S, and let v be the column of index j of the matrix /5(k). The only 
non-zero entry of v, which is 1, is located at the index I such that 1 = kj, i.e. [n/l] = [n/jk]. 
Therefore, the column Tv is the column of index I of T, which is composed of I's for all indices i 
such that i <l = [n/l], and of O's below. Moreover, the integers / and I are in symmetric positions 
in the list S (see Proposition 4). From this we deduce that 

{Tp{k)). . = 1 ^ {Tv)i = 1 ^ i< [n/l] ^ i < [n/jk] 
<^ i < n/jk <^ ij < n/k <^ ij < [n/k] <^ k < [n/ij] 

Proposition 20. For all a G .A, the matrix Tp{a) is symmetric. 

Proof. From Lemma 19, the matrices Tp(k) are symmetric. Moreover, these matrices form a 
basis of A.. Therefore, by linearity, the matrix Tp{a) is symmetric for every a G -4. 

Proposition 21. If we introduce the notations U = Tp{u) and M = Tp{yb), then the matrices U 
and M. are symmetric and satisfy the relation 

M = TU'^T. 

Proof. From Proposition 13 and item 3 of Proposition 9, we have /x = u"-*^, from which it follows 
that 

M = Tp{n) = Tp{u)-^ = T (U-^T) . 
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2.5. The matrix Ai consists of values of the Mertens function 

Proposition 21 establishes a relation between the matrices U and Ai. In essence, we can 
compute M by inverting U and multiplying the result on both sides by T. The following is 
another relation between U and M. , involving the Mertens function. 

Proposition 22. For every n eW , the matrices U and M can be computed by the formulas 

U = {[n/ij]),j^s ''^d M = {M {[n/ij]))^^.^s ' 

where M denotes the Mertens function. In other words, M can be computed by applying the 
Mertens function term by term to the matrix U. 

Proof. We noted in Corrolary 14 that u = Ylkesi^ ~ k~)\^. Therefore, by linearity, 

U = Ekesik-k-)Tpik), 

and from Lemma 19, we deduce that 

^hj= XI ik-k-) = [n/ij]. 

keS, k<[n/ij] 

This completes the proof concerning u. 

Similarly, from Corollary 14, we have /x = J2keS (-^(^) ~ M{k~)) k. Hence, 

M = Z,^smk)-M{k-))Tp{k), 

Mij = Y.kesmk)-M{k-)){Tp{k)\., 

Mij = Ekes, k<[n/^,] mk) - M{k-)) , 

Mi J = M{[n/ij\). 

Example 23. For n = 16, the matrices U and M are 



u 


1 


2 


3 


4 


5 


8 


16 


M 


1 


2 


3 


4 


5 8 


16 


1 


16 


8 


5 


4 


3 


2 


1 


1 


-1 


-2 


-2 


-1 


-1 


1 


2 


8 


4 


2 


2 


1 


1 




2 


-2 


-1 








1 1 




3 


5 


2 


1 


1 


1 






3 


-2 





1 


1 


1 




4 


4 


2 


1 


1 








4 


-1 





1 


1 






5 


3 


1 


1 










5 


-1 


1 


1 








8 


2 


1 












8 





1 










16 


1 














16 


1 













Throughout this section, the matrices M. , T and U were not indexed by the integer n, although 
these matrices depended on n. Henceforth, we will use the notations r„, Un, Mn instead of T, U, 
M, to express the dependence of these matrices on n. 

Theorem 24. The Riemann hypothesis is true if 

\\Mn\\ = 0(n^/^+^), for every e > 0. 

Proof. On the one hand, we derive from Proposition 22 that M{n) is the (1, l)-entry of Mn, and 
we know that for every square matrix A, one has max|^jj| < ||A|| (see [5] p. 57). Therefore, we 
have |M(n)| < ||A^n||- On the other hand, Littlewood [11] proved that the Riemann hypothesis 
is equivalent to the estimate M{n) = 0(n^/^"'"'^), for every e > 0. 



8 



Remark 25. Our hope is that a good estimate of the spectral radius of A4„, i.e. its largest 
eigenvalue, could eventually be obtained by investigating the smallest eigenvalue of the inverse 
matrix T~^Z^„T~^, whose construction is quite simple. 

In the next section, we will look at the quantity ||A^n|| experimentally, as n varies through a 
certain range of integers. 



3. Numerical experiments 

This section presents the results of some experimental computations concerning the growth of 
the sequence ||A^n|| a-s n tends to infinity. 

3.1. Regularity of the sequence \\Airi\\ 

Figures 1 and 2 display the sequences M{n)lyjn and ||A^n||/\/^ respectively, for n running 
from 10^ to 10^, with a step of 10^. Figure 2bis shows the same data as in Figure 2, but displayed 
in a window of smaller height. 



Figure 1 




2 4 6 8 10 



" x10' 
Figure 2bis 

1 .305 I \ \ \ \ 1 




We observe that the growth of ||7W„||/-y/n is quite regular, in contrast to the chaotic behavior 
of M{n)l \/n. Not only is the behavior of ||A^„||/-y/n more regular, but Figures 2 and 2bis also 
show that the range in which the sequence takes its values is much narrower as n 

increases, as compared to the case of M{n) / ^/n. Another important observation is that the 
growth of ||A4„||/y^ seems to be relatively slow. We now look more closely at this growth. 
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3.2. Experimental convergence of the sequence — ^[^^"^"^^^ towards 1/2 

log n 

With regard to Theorem 24 we now turn our attention to the sequence Wn = — ^^H"^"!!-^ — 1/2. 

log n 

Figures 3 to 6 display the sequence Wn for n taking all the integer values in four intervals centered 
on the values n2 = 200^, ns = 300^, = 400^, and n5 = 500^. 



0.0246 



0.0246 



0.0245 



3.95 



Figure 3 




Figure 5 



X 10 




0.023 



Figure 4 




8.95 



Figure 6 



2.49 



X 10 



2.5 
n 



9.05 
X 10" 




2.51 
X 10^ 



In these figures, the points (n, Wn) are linked by a blue line, but when n is of the form k"^ or 
fe^ + k, we plot (n, Wn) as a red diamond. The diamonds are linked by a solid line for greater 
legibility. We distinguish the two cases because the size of the matrix increases by one pre- 
cisely when n increases from /c^ — 1 to k"^ or from A;^ + /c — 1 to /c^ + /c. On the figures, this results 
in small upward jumps in the values of Wn- Some structure can be observed in the variations of 
Wn between every two successive integers n of the form k"^ or k"^ + k. Inside such an interval, the 
values of Wn seems to follow a random walk of moderate amplitude with a dominant decreasing 
trend, followed by jumping when n reaches k"^ or k"^ + k. These observations suggest that the 
overall behavior of the sequence Wn is best described when the values of n are restricted to the 
forms k"^ or k"^ + k. 



Figure 7 displays the sequence Wn for n running from 10'^ to 10^, with the values of n restricted 
to the forms /c^ or k'^ + k] Figure 8 shows the same data displayed in loglog axes, i.e. \og{wn) 
plotted against log(n). 
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Figure 7 Figure 8 

0.034 
0.032 



0.03 



0.028 



0.026 



0.024 



0.022 



0.02 
0.018 

5 10 6 8 10 12 14 

n xlO^ 

We observe in Figure 7 that Wn is roughly decreasing and remains positive (the positivity 
results from the fact that ||A4„||/-y/n > 1 within the range considered), so we may expect that 
Wn converges to some limit, possibly 0. The possible convergence of Wn towards is not very 
apparent in Figure 7, but it is more striking in Figure 8. If, as the graph suggests, this trend were 
to be confirmed as n increases indefinitely, then the following conjecture would be true : 

Conjecture 26. We have the estimate : 

\\Mn\\ = 0(n^/2+'), for every e > 0. 

4. Conclusion 

We have built a sequence of symmetric matrices Mn satisfying |M(n)| < for all pos- 

itive integers n, where M denotes the Mertens function. Based on numerical evidence, we have 
conjectured that ||A1n|| = 0{n^^'^~^^), for every e > 0, a statement which implies the Riemann 
hypothesis. It may be noted that in no part of this study have we made use of complex variable 
methods. Finally, the symmetry of the matrices A4n suggests that spectral methods in matrix 
analysis could play a more significant role in the search for a solution to the Riemann hypothesis. 
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